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ABSTRACT 

We present results from integral field optical spectroscopy with the Potsdam Multi- 
Aperture Spectrograph of the Herbig-Haro (HH) object HH 204, with a spatial sam- 
pling of 1 x 1 arcsec 2 . We have obtained maps of different emission lines, physical 
conditions and ionic abundances from collisionally excited lines. The ionization struc- 
ture of the object indicates that the head of the bow shock is optically thick and has 
developed a trapped ionization front. The density at the head is at least five times 
larger than in the background ionized gas. We discover a narrow arc of high T e ([N n]) 
values delineating the southeast edge of the head. The temperature in this zone is 
about 1,000 K higher than in the rest of the field and should correspond to a shock- 
heated zone at the leading working surface of the gas flow. This is the first time this 
kind of feature is observed in a photoionized HH object. We find that the + and O 
abundance maps show anomalous values at separate areas of the bow shock probably 
due to: a) overestimation of the collisional de-excitation effects of the [O n] lines in the 
compressed gas at the head of the bow shock, and b) the use of a too high T e ([N II]) 
at the area of the leading working surface of the flow. 

Key words: ISM: abundances - ISM: Herbig-Haro object - ISM: individual: Orion 
Nebula - ISM: individual: HH 204 



1 INTRODUCTION 

Herbig-Haro (HH) objects were originally discovered as 
knots of optical emi ssion, found i n regions of low- 
mass star formation jHerbid Il95ll : iHarcT ll952T ). How- 
ever, they are now observed in all spectral domains 
and ar e also associated with high-mass star-forming re- 
gions l|Marti. Rodriguez fc Reipur th 1993). They are of- 
ten found in groups, arranged in linear or quasi-linear 
structures, frequently with aligned proper motion vectors 
jHerbig fe Jones! Il98ll: iReipurth. Raga fe Heathcotel Il99 j ; 
Eislof fel. Miindt fe Bohmlll994l 1 and are a manifestation of 
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the collimated supersonic ejection of material from young 
stars or their close circumstellar environment. HH objects 
show an emission line spectrum characteristic of shock- 
waves, with velocities typically in the range 50 — 200 
km s" 1 (|Schwartz fe Dopitai ri980') and, in many cases, in- 
dividual knots or groups of knots have shapes that are very 
suggestive of bow shocks. 

There are many HH objects identified in the Orion Neb- 
ula - the nearest high- mass star-forming region- which is the 
most observed and studied Galactic H II region. The most 
prominent high-velocity feature in the nebula is the Becklin- 
Neugebauer/Kleinmann-Low (BN/KL) complex, which con- 
tains several HH objects. In addition, there are other impor- 
tant high- velocity flows that do not belong to the BN/KL 
complex, as is the case of HH 202, 203 and 204. The origin of 
these flows has been associated wit h infrared sources embed- 
ded within the Orion- South (see lO'Dell fc Hennevl l2008al . 
and references therein). The hottest star in the Trapezium 
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Figure 1. Left-hand s ide: HST image taken with the Wide-Field Planetary Camera 2 (WFPC2) in the narrow filter F656N (Hoj+[N II]) 
jQ'Dell fc Wonelll99&t) . The large white square corresponds to the usable PMAS FOV covering the Herbig Haro object HH 204. Right- 
hand side: PMAS Ha emission map with HST contours overplottcd. 



cluster, 9 1 Ori C (07 V. ISimon-Dfaz et al]l2006h . dominates 
the photoionization of the surrounding gas and produces a 
thin blister of photoionized material. In general, the mixture 
of photoionized and shocked gas makes the interpretation of 
the spectra of HH objects within H II regions more difficult; 
for example, emission lines formed as a shock moves through 
ionized gas clos ely resemble those from higher velocity shock 
in neutral gas (|Cox fc Ravmondlll985l ). Understanding line 
ratios in H II regions also requires distinguishing shocked 
from photoionized gas before we determine canonical elec- 
tron temperatures and densities and ionic abundances of 
such regions. 

HH 204 -as well as th e neighbouring HH 203- is located 
close to 8 2 Ori A (09 V. ISimon-Dfaz et alJkood ) and has 
an optical angula r size of the order of tens of arcseconds. It 
was discovered bvlMunch fc Wilson! l| 19621 ') and identified as 
HH object by ICanto et all |l980l ). It has been the subject 
of numerous subsequ ent studies of their ionization st ructure 
dO'Dell et alJll997bT). electron den sity (|Walshlll982T ). radial 
velocity dO'Dell. Wen fc Hulll993r) and tangential motions 
JCudworth fc Stone! 1*9771; iHull 19961 ; Iboi. O'Dell fc Hartiganl 
l2002MRosado et all 1)2002 )7 in a kinematic study based on 
Fabry-Perot data, suggested that it seems to be part of a 
large structure or lobe. Also, these authors showed that the 
[O ill] emission does not come from the head of the bow 
shock, but instead occurs primarily in a diffuse region along 
the southwestern portion of the object. lO'Dell et al.l (|l997a ) 
attribute the distribution of [O ill] to differences in the way 
in which 8 1 Ori C illuminates this region of the nebula. The 
ultraviolet light from this star penetrates into the back side 
of the shock which moves toward the observer ijO'Dell et al.l 
1 1997b! ) 

HH 204 is a highly filamentary bright knot at the end 



of a rather s ymmet ric cone that resembles a bow shock. 
lO'Dell et all (l997al ) indicate that the excitation in this ob- 
ject appears to be controlled by a mixture of local shocks 
and external ionization from 6 1 Ori C. It presents a pro- 
nounced asymmetry in brightness between the bow shock 
side near 9 2 Ori A and the side located away from t his star, 
which is better appreciated in [N n] than in Ha. iHennevI 
119961) has proposed that a transverse density gradient in 
the ambient medium where a bow shock propagates could 
lead to an asymmet ry in brightness of t he bow shock. This 
idea is supported bv lRosado et alj (|200ll ) , who found a slight 
velocity gradient running perpendicular to the axis of HH 
204 in the sense that the fainter (and less de n se) reg ions 
show larger velocities. Originally, lO'Dell et al.l (|l997bl) in- 
terpreted that the shock of HH 204 was forme d when the je t 
hits neutral material in the foreground Veil (jO'Delll 1200 it ) . 
but more accurat e knowledge of the position of Orion's Veil 
l|Abel et al.ll2004l ) and an im proved determination of the tra - 
jectory of the shock allowed iDoi. O'Dell fc Hartiganl (|2004l ) 
to establish that HH 204 arises where the jet strikes denser 
nebular material located where the main ionization front of 
the Orion Nebula tilts up. This tilt is what gives rise to the 
Bright Bar feature that runs from northeast to southwest, 
passing slightly at the northwest of 9 2 Ori A. 

In previous papers, our group has studied the spatial 
distribution and properties of the nebular gas at small spa- 
tial scales in the Orion Nebula. Using long-slit spectroscopy 
at spatial scales of 1.2 arcsec and crossing different morpho- 
logica l structures, iMesa-Delgado. Esteban fc Garcfa-Roiasl 
(2008) found spikes in the distribution of the electron 
temperature and density which are related to the posi- 
tion of proplyds and HH objects. In addition, we have 
published two additional studies using integral-field spec- 
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troscopy at spatial scales of 1 x 1 arcsec 2 , focused on 
the analysis of certa i n inter esting morphological structures. 
iMesa-Delgado et al.l (l2009bh st udied the prominent HH 202 
and Mesa-Delgado et al.l |20l"lT ) analysed areas of the Bright 
Bar and the northeast of the Orion South Cloud. They have 
mapped the distribution of the emission line fluxes and the 
physical and chemical properties of the nebular gas in the 
selected areas. 

The main goal of this paper is to use integral field spec- 
troscopy at spatial scales of 1 arcsec in order to explore the 
ionization structure, physical conditions and chemical abun- 
dances of one of the most remarkable HH object of the Orion 
Nebula: HH 204. 

In Section 2 we describe the observations obtained and 
the reduction procedure. In Section 3 we describe the dif- 
ferential atmospheric refraction correction, the emission-line 
measurements and the reddening correction of the spectra 
as well as the determination of the physical conditions and 
chemical abundances. In Section 4 we present the spatial 
distribution of emission line ratios, physical conditions and 
chemical abundances. In Section 5 we discuss the trapped 
ionization front and the high-temperature arc. Finally, in 
Section 6 we summarize our main conclusions. 



1.0 



0.5 



0.0 



-0.5 



-1.0 



H I lines n~ 
He I lines a" 



4000 



5000 6000 
Wavelength (A) 



7000 



Figure 2. Comparison between the differential atmospheric re- 
fraction (DAR) values computed from Filippcnko's model (solid 
lines) and the declination shifts obtained from our data using the 
CORREL.OPTIMIZE routine avaible in the IDL Astronomy User's 
Library. The squares represent the shifts of the H I lines with re- 
spect to H/3 (filled) and the triangles the He I lines with respect 
to He I 4921 A (filled). The gap in the continuous line shows the 
part of the optical spectrum do not covered with the blue and red 
ranges of our observational setup. The symbols in the red range 
have been shifted —0.20 arcseconds (see text). 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were made on 2008 December 20 at 
the 3.5m Calar Alto telescope (Almen'a, Spai n) with the 
Postd am Multi- Aperture Spectrometer (PMAS. lRoth et al.l 
120051 ) in service mode. The standard lens array integral field 
unit (IFU) was used with a sampling of 1" and 16" x 16" 
field of view (FOV). The blue range (3500 to 5100 A) and red 
range (5700 to 7200 A) were covered with the V600 grating 
using two grating rotator angles: —72 and —68, respectively, 
with an effective spectral resolution of 3.5 A. The blue and 
red spectra have a total integration time of 2400 and 300s, 
respectively. Additional short exposures of 30s were taken 
in order to avoid saturation of the brightest emission lines. 
Calibration images were obtained during the night: contin- 
uum lamp needed to extract the 256 apertures, arc lamps for 
the wavelength calibration and spectra of the spectrophoto- 
metric stan dard stars Feige 34, Feige 110 and G 191-B2B 
|Oke]|l99rJ ) for flux calibration. The error of the absolute 
flux calibration of the spectra is of the order of 5%. The 
observing night had an average seeing of about 1.5". 

The data were reduced using the iraf reduction pack- 
age SPECRED. After bias subtraction, spectra were traced 
on the continuum lamp exposure obtained before and after 
each science exposure, and wavelength calibrated using a 
Hg-Ne arc lamp. The dome flats were used to determine the 
response of the instrument for each fiber and wavelength. 
Finally, for the standard stars we have co-added the spectra 
of the central fibers and compared them with the tabulated 
one-dimensional spectra. 

On the left-hand side of Fig. [1] we show a high resolu- 
tion Hubble Space Telescope (HST) image taken with the 
Wide-Field Planetary Camera 2 (WFPC2) in the narrow- 
band filter F656N (Ha + continuum) retrieved from its pub- 
lic archive. Two images were combined using the task crrej 
of the hst_calib.wfpc package of iraf to reject cosmic rays. 
To locate precisely the PMAS FOV on the HST image, we 



have used the star marked with a small white box in Fig. \T\ 
(left), which is visible in the spaxel (—2,6) in all the con- 
tinuum maps of the PMAS field. Therefore, the accuracy in 
the location of the usable PMAS FOV (the large white box 
on the HST image) is better than 1". With the //"ST image, 
we computed high-resolution contours and overplotted them 
on all maps. Note that the resolution of the HST image is 
roughly one order of magnitude better than the PMAS one 
(see right-hand side of Fig. [TJ . 

3 ANALYSIS 

3.1 DAR correction 

The effect of the differential atmospheric refraction (DAR) 
can be noticed in the data. We have used the monochro- 
matic images of H I and He I emission lines at different 
wavelengths -using H/3 and He I 4921 A as references- to 
correct for this effect. With the CORREL_optimize routine 
avalaible in the IDL Astronomy User's Library we computed 
the optimal pixel offset of the monochromatic image of each 
selected line relative to its reference by maximizing the cor- 
relation function. In Figure[2]we represent the displacements 
in declination due to t he DAR compare d with the results 
of Filippenko's model (|Filippenkd Il982h and with respect 
to H/3 and He I 4921 A. As the optical spectra analyzed 
are separated in two spectral ranges, there is an additional 
fixed instrumental shift between the lines of the red and blue 
ranges. This shift has been double checked using the pair of 
images Ha— H/3 and He I 4921 - 6678 A, yielding displace- 
ments of (0.10, —0.20) and (0.10, —0.25) pixels, respectively. 
We chose the first pair of values because it is computed with 
much brighter lines. Therefore, the displacements of lines in 
the red range of the spectrum have been shifted (0.10, —0.20) 
arcsecs. Figure [2]only presents the behaviour of DAR in the 
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Figure 3. Spatial distribution map of the extinction coefficient, 
c(H/9). The contours plotted are those from the HST F656N im- 
ages. The size of the field corresponds to the maximal coincident 
FOV after correcting the spectra for differential atmospheric re- 
fraction (see text). 



declination axis because it has larger values than in right 
ascension due to the paralactic angle is better aligned with 
the declination axis. Attending to the good agreement of 
the model with the data values, we corrected from DAR the 
whole set of data cube applying the bilinear routine of idl 
with a global and independent offset to each axis, in order 
to maximize the coincident FOV (15 x 15 arcsec 2 ). 



3.3 Reddening coefficient 

The reddening coefficient, c(H/3), has been obtained by fit- 
ting the observed H/y/H/3 and H^/H/3 ratios t o the the- 
oretical ones predicted by IStorev fc Hummer] l|l995h for 
n e = 1000 cm' 3 and T e = 10000 K. We have used 
the redd ening function, /(A), normalized to H/3 deter- 
mined by iBlagrave et all (120071 ) for the Orion Nebula. The 
use of this extinction law instead of the classical one of 
ICostero fc Peimbertl ||l970l ) produces c(H/3) values about 
0.1 dex higher and dereddened fluxes with respect to H/3 
about 3% lower for lines in the range 5000 to 7500 A and 
4% higher for wavelengths below 5000 A. The final adopted 
c(H/3) value for each spaxel is an average of the individ- 
ual values derived from each Balmer line ratio weighted by 
their corresponding uncertainties. The typical error is about 
0.1 dex for each spaxel. The resulting extinction map is 
shown in Figure [3] where the extinction coefficient varies ap- 
proximately from 0.1 to 0.6 dex. The mean value in the FOV 
is 0.23 dex. We have com pared our c(H/3) value with t he ex- 
tinction map obtained bv lo'Dell fc Yusef-Zadeh] (|2000l ) from 
the Ha/H/3 line ratio using calibrated HST images, where 
HH 204 is located between the contour lines of 0.2 and 0.4 
dex on that map. Therefore, we find that our mean value 
of the reddening coefficient is consistent with the previous 
determinations taking into account the errors and the dif- 
ferent extinction laws used. However, we also find some re- 
gio ns with higher values o f exti nction than those reported 
by lO'Dell fc Yusef-Zadehl (|2000l ). ft must to be said that, 
due to its closeness to 8 2 Ori A, the H5 and H7 lines of the 
spectra are somewhat affected by underlying absorption due 
to the stellar continuum reflected by dust. This may affect 
the measurement of the intensity of Balmer lines, especially 
the bluest ones, producing a slightly higher c(H/3), and this 
may be the reason of such discrepancy. However, the results 
of our paper are not particularly sensitive to this issue. 



3.2 Line intensity measurements 

The emission lines considered in our analysis are the follow- 
ing: hydrogen Balmer lines, from Ha to H12, which are used 
to verify the DAR correction and to compute the redden- 
ing coefficient; collisionally excited lines (CELs) of various 
species, which are used to compute physical conditions and 
chemical abundances. 

Line fluxes were measured applying a single or 
a multiple Gaussian profile fit procedure between two 
given limits and over the local continuum. All these 
measurements were made with the SPLOT routine of 
1 RAF and using our own scripts to automatize the pro- 
cess. The associated errors were determin ed following 
iMesa-Delgado. Esteban fc Garcia-Roiasl (|2008l ). In order to 
avoid spurious weak line measurements, we imposed three 
criteria to discriminate between real features and noise: 1) 
Line intensity peak over 2.5 times the sigma of the contin- 
uum; 2) FWHM(H i)/1.5 < FWHM(A) < 1.5xFWHM(H 1); 
and 3) F(A) > 0.0001 x F(H/3). 

All line fluxes of a given spectrum have been normalized to 
H/3 and Ha for the blue and red range, respectively. To pro- 
duce a final homogeneous set of line ratios, all of them were 
re-scaled to H/3. The re-scaling factor used in the red spectra 
was the theoretical Ha/H/3 ratio for the physical conditions 
of T e = 10000 K and n e = 1000 cm" 3 . 



3.4 Physical conditions 

Nebular electron densities, n e , and temperatures, T e , have 
been derived from the usual CEL diagnostic ratios - 
[S 11] 6717/6731 for n e , and [N 11] (6548+6583)/5755, [S 11] 
(6716+6731)/(4068+4076) and [O III] (4959+5007)/4363 
for T e - and using the temden task of the neb- 
ula r package of iraf with updated atomic data 
(see [G arcia- Roias. Pefia fc Peimbertl 1 2009 ) . Following the 



same procedure as Mesa-Delgado. Esteban fc Garcia-Roiasl 
(2008), we have assumed an initial T e — 10000 K to de- 
rive a first approximation of n e ([S 11]), then we calculate 
T e ([N 11]), and iterate until convergence. We have not cor- 
rected the observed intensity of [N 11] 5755 A for contribu- 
tion by recombination when determining T e ([N 11]) because 
this is expected to b e rather small in the Orion Nebula (e.g. 
lEsteban et~al"1l2004l ). We did not use T e ([S 11]) for the itera- 
tion procedure because the auroral lines of [S 11] are rather 
weak outside HH 204 and the diagnostic line ratio involves 
lines located in the blue and red spectral ranges. Neverthe- 
less, using both [S 11] diagnostic ratios together the density 
values are in good agreement within the errors with those 
determined with T e ([N 11]). On the other hand, because of 
the relatively low spectral resolution of the observations, 
[O in] 4363 A is somewhat blended with the Hg l+[Fe 11] 
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4358 A lines. The latter one becomes very strong inside the 
HH object. Only in this case, we performed manually the 
Gaussian fit of those three lines in order to obtain a good 
determination of T e ([0 ill]). The automatic procedure gave 
some inconsistent values at the northeast corner of the FOV 
because of the faintness of the [O ill] 4363 A line in that 
area. The typical uncertainties in the density maps range 
from 200 to 2000 cm -3 . The errors for electron tempera- 
tures range from 150 to 350 K for T e ([N II]), from 100 to 
400 K for T e ([0 ill]) and from 100 to 500 K for T e ([S n]). 

3.5 Chemical abundances 

The iraf package nebular has been used to derive ionic 
abundances of S + , N + , + and 2+ from the intensity of 
CELs. We have assumed no temperature fluctuations in the 
ionized gas (t 2 — 0) and a two-zone scheme: T e ([N n]) is 
used to derive the S + , N + and + -the low ionization po- 
tential ions- and T e ([0 ill]) for the 2+ abundances. The 
uncertainties in ionic abundances have been calculated as 
a quadratic sum of the independent contributions of errors 
in n e and T e . For oxygen abundances, the uncertainties are 
0.04-0.08 dex for 0+, 0.02-0.05 dex for 2+ and 0.02-0.06 
dex for O. 



4 RESULTS: SPATIAL DISTRIBUTION MAPS 

In this section we analyse the spatial distributions of several 
emission line ratios, electron densities and temperatures as 
well as several representative ionic and total abundances. 

4.1 Ionization structure 

In Fig. 01 we show the line intensity ratio maps of [O i] 6300 
A, [S n] 6731 A, [O n] 3727 A, [N n] 6583 A, [Fe in] 4881 
A, [S in] 6312 A, He I 6678 A, [O ill] 5007 A and [Ne ill] 
3868 A relative to K/3. The minimum and maximum values 
of the color bars of the different ratios have been rescaled to 
show the same dynamical range in all the maps. The dynam- 
ical range of reference adopted is that of the [O i]/H/3 ratio, 
which shows the largest variation. The selected lines cover 
all possible ionization conditions of the gas. The maps are 
sorted from lower (top-left) to higher (bottom-right) ioniza- 
tion potential. These line ratios are directly related to the 
local ionization degree of the gas and do not depend on den- 
sity in a direct manner. 

[O i] 6300 A is a strong l ine in shocks, but is rathe r weak 
in photoionized regions (e.g. iMasciadri fe Ragall200ll ) . How- 
ever, this line can also be a tracer of the presence of trapped 
ionizatio n fronts (as in the case of HH 202, in the Orion 
Nebula: iMesa-Delgado et al.l [2009al lbl^ ■ In Fig. QJi, we can 
see that [O i]/H/3 is much higher in the southeast quadrant 
of the bow shock and very low along the bow shock wings. 
This indicates that the wings of the flow are completely pho- 
toionized. In order to disentangle the excitation mechanism 
of the gas at the bow shock, we have represented its [O l]/Ha 
and [O l]/[Q hi] ratios on a diagnos tic diagram adapted from 
iRiera. Phillips fe Mampasd (|l990l ) that separates shock ex- 
citation from photoionization, finding that it is esentially 
photoionized, and that no substantial contribution of shock 



excitation is needed to reproduce its spectrum. In their mod- 
els of optically thick HH objects inside photoionized regions, 
IMasciadri fe Ragal (|200lh find that trapped ionization fronts 
restricted at the head of the bow shock show an [O i]/H/3 ra- 
tio concentrated to the side directed away from the ionizing 
photon source, exactly what we see in Fig. The [S il] line 
emission is also a tracer of ionization fronts. In Fig. |4p, wc 
can see that [S ll]/H/3 follows rather closely the distribution 
of[Oi]/H/3. 

As we can see in Fig. [3t and [IJi, the distribution of 
[O ll]/H/3 shows an almost inverse behaviour with respect to 
[O mJ/H/? presenting its lowest values at the western part 
of the FOV, exactly where the [O mJ/H/? ratio is higher. 
However, [O m]/H/3 reaches its minimum at the head of the 
bow shock while [O ll]/H/3 does not present peak values at 
those spaxels. This indicates that the spatial distribution of 
[O II] 3727 A and H/3 is rather similar at the bow shock. 

The [N ll]/H/? ratio (Fig. 2J1) reaches its maximum at 
the bow shock, and delineates the wings of the flow struc- 
ture. It can also be seen that the area closer to the star 
6 2 O ri A is brighter than the rest. This result is in agreement 
with lRosado et all l|200lf ), who reported that the asymmetry 
in the brightness distribution of HH 204 is better appreci- 
ated in [N n] than in Ha. 

HH objects are also characterized by their strong emis- 
sion in [Fe ill] lines. In Fig.0^, we can see the spatial distri- 
bution of the [Fe iii]/H/3 ratio where the maximum values 
seem to be concentrated just behind the areas of the max- 
ima of the [O l]/H/3 and [S ii]/H/3 ratios. The high [Fe ill] 
emission may be related to destruction of dust grains due 
to the passage of shock w aves as in the case of HH 202 
jMesa-Delgado et al.|[2009aH bl). 

The spatial distributions of [S iii]/H,5 and He l/H/3 
are rather similar and smooth (Fig. |4f, g)- Higher values 
tend to be concentrated toward the west of the FOV. This 
concentration is more evident in the higher ionized species 
[O m]/H/3 and [Ne m]/H/3 (Fig. gh, i) due to the contrast 
between the minimum and maximum line ratios is more pro- 
nounced. It is important to comment that the high excita- 
tion lines of [O ill] and [Ne ill] do not form at the head of 
the bow shock like they usually do in other HH objects pro- 
duced by high- velocity flows, but instead occur primarily in 
the ambient diffuse background across the western portion 
of the FOV. This distribution is attributed to differences in 
the way in w hich the star 9 1 Ori C illuminates this zone 
of the nebula (|Q'Dell et alJll997bh and the relatively large 
distance between the star and HH 204, which makes the im- 
pinging ionizing flux too low to maintain the presence of a 
substantial fraction of 2+ ions at the bow shock. 

4.2 Physical conditions 

In Fig. [5] we present the physical conditions determined for 
the observed field, where the head of the bow shock is asso- 
ciated with local peaks of electron density. 

Fig- [5^1 shows the T e ([N n]) map. A remarkable feature 
is the narrow high-temperature zone just at the southeast 
border of the head of the bow shock. This structure seems to 
be real, since the temperature determination has only an er- 
ror of about 2%, much lower than the temperature increase. 
Furhtermore, it is a rather narrow feature that can only be 
distinguished when the DAR correction is applied. On the 
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Figure 4. Maps of line intensity ratios relative to H/3 of selected emission lines: a) [O i] 6300 A, b)[S n] 6731 A, c) [O n] 3727 A, d) 
[N II] 6583 A, e) [Fe in] 4658 A, f) [S III] 6312 A, g) He I 6678 A, h) [O in] 5007 A and i) [Ne in] 3868 A. The contours plotted are 
those from the HST F656N images. The color bars have been re-scaled to have the same dynamical range between the maximum and 
minimum values in all maps. 



other hand, this high-temperature zone cannot be produce 
by the slight extinction excess we find in a nearby -but not 
coincident- area of the head. We would need c(H/3) differ- 
ences by about unity to account for a temperature increase 
of 1000 K. We also do the exercise of correcting the whole 
set of data with two constant values of c(H/3), 0.2 and 0.4 
dec, and the high-temperature arc is still clearly present in 
both cases. Then the derived temperature enhancement is 
not an artifact of the high local extinction. The nature of this 
structure will be discussed later. Fig.[5p shows the n e ([S n]) 
map derived in combination with T e ([N n]) until conver- 
gence. Assuming this initial n e ([S n]) map we also derived 
the T e ([S n]) map (see Fig. [5^), which shows an arc that 
delineates a zone of strong increase of temperature at the 
northeast of the head. We found that this structure also ap- 
pears if we calculate the map of the S + abundances derived 
from the auroral [S n] lines when using T e ([N n]) (Fig. [5^) 
and n e ([S n]) (Fig. [SJa). Obviously, that result has no phys- 
ical meaning. The simpler explanation for this fact is that 
we are not using appropriate values of the electron density 
for deriving S + abundances and this affects especially the 
calculations based on nebular [S n] lines, which have lower 



critical densities than the auroral ones. In order to correct 
for this discrepancy in the S + abundance, we determined 
the n e ([S n]) map that -assuming the T e ([N n]) distribution 
of Fig. [5^,- makes the differences between the abundance 
obtained from the nebular and auroral lines lower than the 
estimated maximum error of this abundance, ±0.08 dex. The 
final electron density map is shown in Fig. [5fi, where only 
half of the spaxels show values different from the previous 
map. It must be said that using an iterative procedure in 
which both T e ([N n]) and n e ([S n]) are changed we do not 
find a substantially different final T e ([N n]) map because 
the [N n] lines involved in the temperature indicator have 
critical densities much higher than the n e ([S n]) values we 
are obtaining for HH 204. In Fig. (SJi we can see that this 
new density map fits much better the Ha contours not only 
inside the head of the bow shock, but also at the location of 
other structures visible in the HST images such as the wing 
at the west of the head. We recomputed T e ([S n]) assuming 
the n e ([S n]) map of Fig. [5J; and obtained the distribution 
showed in Fig. [Sf, finding a much smoother map where the 
arc-shaped area of very high temperatures has disappeared. 
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Figure 5. Physical conditions of HH 204 a) T e ([N II]), b) and c) two different estimates of n e ([S II]), d) T e ([0 in]), e) and f) two 
different estimates of T e ([S II]). See text for details. The contours plotted are those from the HST F656N images. In the density map b) 
we mark three different zones: background gas (white box), head of the bow shock (black box) and high-temperatutre arc (grey box), 
whose spaxels have been co-added in order to analyse their spectra and derive the physical conditions at those zones. 



The final n e ([S n]) map of Fig. [5J: is the one we have used 
to derive the chemical abundances. 

In Fig. [5}j we can see that the highest densities in the 
FOV are be tween 4000 and 9000 cm" 3 . Similar values were 
obtained bv lMesa-Delgado. Esteban fc Garcfa-Roiasj (|2008l ) 
from a slit crossing HH 204 and roughly parallel to the di- 
rection toward 8 1 Ori C (see Fig. [7}. The density map shows 
the highest values at the roundish clumpy structure of the 
head of the bow shock while the extended wings show densi- 
ties only slightly higher than those of the surrounding back- 
ground gas. The lowest density values (n e < 2000 cm" 3 ), 
related to the background gas, are located at the southwest 
edge of both FOV. We obtained an additional density deter- 
mination making use of the [Fe ill] lines. We detected at least 
5 [Fe in] lines of the 3F multiplet ([Fe ill] 4607,4658, 4702, 
4734 and 4755 A) and one line of the 2F multiplet ([Fe in] 
4881 A) at those spaxels at the head of the bow shock show- 
ing the maximum values of n e ([S n]). The range of valid- 
ity of this density indicator, n e ([Fe III]), has an upper limit 
well above 10 6 cm -3 . To carry out these calculation we have 
used a 34-level atom using the collision st rengths from Zhang 
Jl996t) and the t r ansiti on probabilities of lQu inct ( 1996]) and 
Ijohansson et al.l (|2000l ). The calculations of n e ([Fe m]) give 
densities varying from 8000 to 15000 cm -3 , but the uncer- 
tainties are rather large due to the faintness of the [Fe in] 
lines. Therefore, the n e ([S n]) values we obtain in Fig. [SJ; 
should be considered a fair approximation to the true elec- 
tron density of the ionized gas at HH 204. It is interesting 



to comment that Solf. Bohm fe Ragal (|l98ct ) found a rather 
similar lack of correlation between auroral and nebular [S n] 
lines in a one-dimensional spectrum across the head of the 
prototypical HH 1. They also found rather high n e ([S n]) 
values -of the order of 3000 cm" 3 - that peak at the head, 
precisely where the intensity ratio between auroral and neb- 
ular [S n] lines is higher. This behavior is also similar to that 
observed in our case and can also be explained by the effect 
of collisional deexcitation of the nebular lines. 

The spatial distribution of T e ([0 ill]) (Fig.[S}l) is rather 
smooth, with higher values located at the northeast edge 
of the FOV. This area shows the faintest [O in] 4363 A 
emission and therefore have the largest errors. In fact, the 
anomalously high T e of some spaxels might indicate that the 
[O ill] 4363 A line has not been properly deblended from 
Hg I+[Fe in] 4358 A lines. There is an apparent trend of 
slightly higher T e ([0 ill]) around the northern half of the 
head of the bow shock. 

In Fig. [5p we mark three selected areas covering dif- 
ferent zones, whose spectra have been extracted in order 
to analyse their physical conditions. These areas are: back- 
ground gas (white box), head of the bow shock (black box) 
and high-temperature arc (grey box). This last zone corre- 
sponds to the area showing the peak of T e ([N n]) located 
at the southeast of the head (see Fig. [5^) . In Table [T] we 
summarize the density and temperature values we obtain for 
each zone. The n e ([S n]) of the three zones are very different, 
1400 cm" 3 for the background, 9000 cm" 3 for the head and 
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Table 1. Physical conditions of the three selected zones (see Fig. [5}}). 



Zone 


c(H/J) 
(dec) 


n e ([S II]) 
(cm -3 ) 


n e ([Fe III]) 
(cm -3 ) 


T=([N II]) 
(K) 


T e ([S II]) 
(K) 


T e ([0 in]) 
(K) 


background 
head 

high-temperature arc 


0.4 ±0.1 
0.3 ±0.1 
0.3 ±0.1 


1400 ± 300 
9000 ± 1800 
3800 ± 550 


2800 ± 1900 
15300 ± 5600 
5000 ± 500 


8800 ± 300 
8700 ± 150 
9600 ± 200 


7700 ± 300 
7800 ± 250 
8500 ± 300 


8500 ± 200 
8700 ± 350 
8800 ± 250 



3800 cm -3 for the high-temperature arc. It is remarkable 
that the electron densities determined with [Fe ill] lines are 
rather consistent with the quoted values of n e ([S n]) within 
the uncertainties and show the same qualitative behaviour. 
On the other hand, the electron temperature of low ion- 
ization species shows similar values in the background and 
the head of the bow shock and higher ones at the high- 
temperature arc. In the case of T e ([N n]) the increase is 
of about 1,000 K and in the case of T e ([S n]) it is some- 
what lower, about 700 K. Finally, T e ([0 m D is practically 
the same in the three areas within the errors. These results 
indicate that there might be a process of extra heating in 
the high-temperature arc area (see Section l5.2p . This pro- 
cess only affects to low-ionization ionic species, because - 
attending to the results of Section and Section |4~3T - high- 
ionization species -such as 2+ - are absent at the head of 
HH 204. 



4.3 Chemical abundances 

In the first row of Fig. [5] we present the spatial distribu- 
tions of the + /H + , 2+ /0 + and O/H ratios derived from 
CELs. The + abundance is strongly concentrated at the 
north of the head (see Fig. [fji). As we can see in Fig. [Ep, 
the 2+ /0 + ratio shows an almost inverse behaviour with 
respect to + abundance, indicating the strong ionization 
stratification across the FOVs. Unexpectedly, we do not find 
a constant value in the spatial distribution of the total O 
abundance, as we can see in Fig. [BJ;. 

The mean value of the O abundance amounts to 
12+log(0/H) = 8.40 ± 0.04 dex, which is slightly lower 
than the typical values of about 8.50 dex obtained 
by previous works in other parts of the Orion Neb- 
ula lEsteban et al.|[i1^.l2004l;lBkgrave, Martin fe Baldwin 
2006; Me sa-Delgado et al.l l2009al ; ISimon-Dfaz fe Stasihska 
20111 ). The maximum values -with abundances between 0.1 
and 0.2 dex higher- are concentrated at the north half of the 
head where the errors are minimal and the minimum values 
-about 0.15 dex lower- are just at the high-temperature 
arc. It is remarkable that the + (Fig. and O (Fig. 0:) 
abundance maps are almost identical, indicating that the 
variations of the total O abundance are produced by the 
structure of the + abundance. A similar result has also 
been obtained in a previous paper of our group devoted to 
the anal ysis of PMAS data of ioniz ation fronts in the Orion 
Nebula l|Mesa-Delgado et al.ll20l"lT ). However, in that paper 
and in contrast to our results for HH 204, we find that the 
0+ and O maps were also similar to that of n e ([S n]). As 
we can see in Fig. [U, the + and O maps reach a minimum 
just at the southeast edge of the head of the bow shock -just 
at the position of the T e ([N n]) peak of Fig. [SJi- and this 



feature does not have any counterpart in neither the initial 
nor the final n e ([S il]) maps (Figs.[5p and c, respectively). 

Considering all the aforementioned results, we propose 
that two different physical effects are affecting the determi- 
nation of the + abundance -and consequently, the total O 
abundance- across th e field of HH 204. Fi r st, fol lowing the 
arguments outlined in lMesa-Delgado et al.l (|201ll ). the most 
likely explanation for the higher + /H + values obtained at 
the head is that the values we obtain for n e ([S n]) are higher 
than the true ones of the + zone. Therefore, we are over- 
estimating the collisional de-excitation of the [O n] 3727 A 
lines. Taking into account t he Q + /H + excess at the h ead 
of HH 204 and the results of iMesa-Delgado et al] (|201ll ) we 
estimate that the true densities at the + zone should be 
about 1000 cm -3 lower than those indicated by n e ([S n]). 
If we correct the density by that factor -which is of the 
order of the quoted errors- we obtain an O abundance sim- 
ilar to that of the background gas around HH 204. Second, 
the effect producing the slightly lower + abundances at 
the high-temperature arc seems to be connected with the 
localized increase of T e ([N n]) seen in Fig. [5^. Assuming 
that the unexpected abundance pattern is the consequence 
of using an electron temperature which is higher than that 
expected considering thermal equilibrium of a static pho- 
toionized nebula, we have estimated by which amount we 
have to correct the T e ([N n]) values of each spaxel of the 
high-temperature arc to reproduce the mean nominal value 
of the total O abundance, 8.40 dex. The result is that, by 
decreasing T e ([N n]) by 1,000 K, we obtain the expected O 
abundance and the arc of high T e ([N il]) disappears. This 
result support the idea that there might be an extra heating 
due to shocks. 

Finally, we present the N+/H+, N+/0+ and N+/S+ ra- 
tios in the bottom row of Fig. \§\ The N + abundance shows 
a rather similar spatial distribution to that of + as it is 
expected from their similar ionization potentials. However, 
the ratio between both ionic abundances shows slight spatial 
variations concentrated at the south edge of the head. These 
variations can be due to the different temperature and den- 
sity dependence of the lines used to derive both ionic abun- 
dances. Finally, the N + /S + ratio delineates the presence of 
a trapped ionization front at the head of the bow shock and 
the strong ionization gradient across the FOV. 



5 DISCUSSION 

5.1 The trapped ionization front 

Most HH objects are located in molecular and neutral envi- 
ronments and their gas is shock-excited. H owever, HH 204 
-like other objects su ch as HH 444 (see iReipurth et al.l 
Il998l ). HH 505 (see iBallv fe Reipurtbl l200ll ). HH ~529 
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Figure 6. Ionic abundance and abundance ratio maps determined from the intensity of collisonally excited lines. From top-left to 
bottom-right: 12 + Iog(0+/H+), log(0 2+ /0+), 12 + log(0/H), 12 + log(N+/H+), log(N+/0+) and log(N+/S+). The contours plotted 
are those from the HST F656N images. 
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Figure 7. Spatial profiles of several parameters along slit position 3 of [Mcsa-Dcleado, Estcb an &; Garcia-Roi crossing HH 204 

on the left panels and a pseudo-slit along the main diagonal of our PMAS maps on the right. Top panels: profiles of n e ([S II]) (solid line) 
and T e ([N II]) (dashed line). Bottom panels: profile of the total oxygen abundance. The origin of abscissae corresponds to the peak value 
ofT e ([Nn]). 



(see Blagrave. Martin & Baldwin! |2006h or HH 202 (see 
Mcsa-Dclgado ct al . 2009a) - is an example of collimated HH 
jet immersed in a photoionized region, where the excita- 
tion of the gas is dominated by the stellar ionizing radia- 
tion. There are rather few theoretical models that describe 
the structure and evolution of photoion i zed HH jets (e.g. 
iRaga et al.ll2000l ; iMasciadri fc Ragall200ll ; iRaga fc Reipurthl 



2004; Raga, Riera fc Gonzalez- Gomez! l20ich . Our unprece- 
dentedly detailed integral field spectroscopical observations 
permit to study not only the ionization structure but also 
the physical conditions of the gas along the head of the bow 
shock, something rather unusual in the literature. Therefore, 
our data are of paramount interest to compare with the pre- 
dictions of jet models. The most suitable models for the par- 
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ticular conditions of HH 204 are those bv lMasciadri fe Ragal 
who consider a source of radiation with a flux of ion- 
izing photon comparable to that expected for 8 1 Ori C and 
located at different distances from the jet. Especially inter- 
esting is their model 2, in which the linear distance between 
the ionizing source and the jet is of the order of that of 
HH 204 with respect to 8 1 Ori C. An apparent drawback o f 
our comparison with the models of lMasciadri fc Ragal (120011 ) 
is that they consider that the direction toward the ionizing 
star is perpendicular to the outflow axis, whilst in our case 
the jet is oriented very close to the direction of 9 Ori C. In 
any case, this difference does not affect the results qualita- 
tively; 

iMasciadri fc Ragal |200ll ) find that their model 2 pro- 
duces a trapped ionization front at the head of the bow 
shock while the wings are completely ionized. This is pre- 
cisely what we observe in our Fig. [4^,, where the [O i] 6300 
A emission is absent in the wings but strong at the head. 
The [O i] emission is concentrated to the south east side of 
the head, roughly directed away from 8 1 Ori C. Evidences of 
a trappe d ionization front were also obtained in the case of 
HH 202 |Mesa-Delgado et al.ll2009al lbT). located at the north 
west of e 1 Ori C. In HH 202, the stronger [O i] 6300 A emis- 
sion is also located in the direction away from the Trapezium 
stars. 

Due to the fact of having a trappe d ionization front and 
following iMesa-Delgado et all l|2009bh . we can estimate the 
width of the ionized slab of HH 204 and its physical separa- 
tion with respect to the main ionization source, 8 1 Ori C. In 
our case, the ionized slab was computed from the maximum 
emission profiles of O and He + with an angular distance 
on the plane of sky of about 2.9 ± 0.5 arcsec. Using t he dis - 
tance to the Orion Nebula obtained bv lMenten et all l|2007f) . 
d = 414 ± 7 pc, and the inclination angle of HH 204 with re- 
spect t o the plane of the sky calculated bv lO'Dell fc Hennevl 
^OOSal ). 9 = 26°, we estimate (6.1 ± 1.2) x 10" 3 pc for the 
width of the ionized slab. To trap the ionization front the 
incident Lyman continuum flux must be balanced by the 
recombination in the ionized slab, resulting in a physical 
separation of D = 0.37 ± 0.10 pc for the n e ([S n]) of the 
head (see Table [T]). This result suggests that HH 204 is em- 
bedded within the body of the Orion Nebula and, therefore, 
discards the origin of the ionized front as result of the inter- 
action _of_the_gas flow_ with the veil, which is between 1 and 
3 pc ( Abel et al. 1 120041 ). If we consider the same ion (from 
O to 0+) to compute the ionized slab ((2.9 ± 1.2) x 10~ 3 
pc), the physical separation is about D — 0.53 ± 0.17. 

5.2 The high-temperature arc at HH 204 

In the models of photoionized HH jets of iRaga fc Reipurtbl 
l|2004l ). the temperature along the HH object is the typi- 
cal of a gas in photoionization equilibrium —about 10 4 K— 
but larger at the leading working surface of the bow shock 
because of shock heating. This high-temperature zone is nar- 
row and precedes the area of high-density shocked gas be- 
hind the working surface. The high-temperature arc we see 
in Fig. [5^ may very likely correspond to that predicted struc- 
ture. In fact, the arc delineates the outer edge of the bow 
shock and precedes the high-density compressed zone that 
forms the head of HH 204 as we see in Fig. [5] The up- 
per limit for the width of the arc is 1", which corresponds 



to a linear dimension of ~6xl0 cm at the distance o f 
the Orion Nebula. The models of lMasciadri fc Ragal (|2001T ) 
—more appropriate for the Orion Nebula conditions— pre- 
dict a high-temperature zone at the leading working surface 
with a width of the order of ~4xl0 15 cm, in a good agree- 
ment with our result. 

There is a previous example of a similar structure in the 
literature. Based on spatia lly resolved echelle spectroscopy, 
Solf. Bohm fc Ragal (|l988T l reported a steep peak of the 
T e ([C i]) just at the leading edge of the bow shock of the 
prototypical HH 1, showing values as high as 2 x 10 4 K 
in a region about 1" wide. Because HH 1 belongs to the 
star-forming region Lynds 1641 —located at the Orion Neb- 
ula Cluster— the linear dimensions of this high-temperature 
zone is rather simi l ar to that we estimate for HH 204. 
ISolf. Bohm fc R aga ( 198£j) finds that the spatial changes of 
electron temperature are far much less pronounced for other 
ions, specially in the case of T e ([N n]), which shows hardly 
any variation. An important difference between HH 1 and 
HH 204 is that the gas excitation mechanism in both objects 
are completely different, pure shock excitation in the case of 
HH 1 and photoionization in the case of HH 204. 

The high-temperature arc at the leading edge of 
HH 204 is a real feature and not an artifact be- 
cause it can also be seen in the l ong-sl it data of 
IMesa-Delgado. Esteban fc Garci'a-Roiasl (|2008T l. In Figs. 
[7^ and b we show the spatial profiles of n e ([S n]) 
and T e ([N n]) (Fig. EK) and the O abundance (Fig. 
7]p) of a small section along the s lit p osition 3 of 
Mesa-Delgado. Esteban fc Garci'a-Roiasl (|2008r i crossing the 
head of HH 204. That slit passes through approximately the 
middle of the high-temperature arc. In Fig. [7^,, we can see 
that the T e ([N n]) peak is 2.5" ahead the n e ([S n]) maximum 
and that the T e ([N n]) peak coincides with a minimum of 
the total O abundance, just the effect observed in our PMAS 
data (see Fig.[U). For comparison, in Figs.Et and d we show 
the spatial profiles of n e ([S n]), T e ([N n]) and O abundance 
obtained from a pseudo-slit -1" wide- obtained along the 
main diagonal of our PMAS maps, which is rather close 
to slit position 3 of lMesa-Delgado. Esteban fc Garci'a-Roiasl 
(2008). The similar behaviour of both sets of profiles is re- 
markable. 



6 CONCLUSIONS 

In this paper we analyse integral field optical spectroscopy 
of the area encompassing the head of the bow shock of 
the conspicuous Herbig-Haro object of the Orion Nebula: 
HH 204. We have mapped emission line fluxes and ratios, 
physical conditions, and ionic and total O abundances at 
spatial scales of 1". 

The maps show a clear ionization stratification across 
the object. The presence of strong [O i] emission at the head 
of the bow shock and its absence at the wings indicate that 
those zones are optically thick and have developed a ioniza- 
tion front. We find evidences that HH 204 is formed within 
the ionized gas of the nebula and not the Veil. The head of 
the bow shock shows electron densities at least five times 
larger than the background ionized gas. We find that the 
density indicator based on [S n] is not providing the cor- 
rect n e ([S n]) values because the nebular and auroral lines 
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of [S n] give different ionic abundances. This is because the 
densities at the head are larger than the critical density. We 
recompute the n e ([S n]) map solving this inconsistency and 
use it in the subsequent calculations. 

We discover a narrow arc of high T e ([N n]) values just 
at the southeast edge of the head. The temperature in this 
zone reaches values about 1,000 K higher than in the rest of 
the object. This arc may correspond to st ructures predicted 
by the m odels of photoionized HH jets of lRaga fe Reipurthl 
(|2004h or lMasciadri fe Ragal l|200ll ). which are produced by 
shock heating at the leading working surface of the gas flow. 
The + and O abundances are somewhat underestimated 
-between 0.1 and 0.2 dex- at this arc due to the higher 
temperature. We find also that the + and O abundances 
are both about 0.1 de x larger at the northern half of the head 
of HH 204. Following [Mesa-Delgado et al.l (|201ll ). the most 
likely explanation is that the appropriate density for the + 
ion should be about 1000 cm -3 lower than the nominal value 
of n e ([S n]) assumed for those spaxels. The use of a wrong 
higher density produces an overestimate of the collisional 
de-excitation of the [O n] 3727 A line. 

The results of this paper demonstrate that the com- 
pression and heating of the gas due to high-velocity flows 
affect the chemical abundance determinations in H II re- 
gions. When densities reach values of the order of 10 3 — 10 4 
cm -3 -as in the c ompressed heads of HH ob jects or at ion- 
ization fronts fsee iMesa-Delgado et al]|201ll )- the + and 
O abundances determined from the intensity of collisional 
excited lines may be overestimated. Also, the presence of 
localized shock heating may affect the thermal balance of 
the gas producing artificial lower abundances. Fortunately, 
these disturbing effects can be analyzed in some detail in the 
Orion Nebula because of the higher spatial resolution attain- 
able for this object due to its closeness, but it is difficult to 
explore in more distant objects and virtually impossible in 
extragalactic H II regions. Similar studies of larger areas of 
the Orion Nebula are necessary in order to: a) quantify how 
these localized effects can affect integrated abundance de- 
terminations when spatially resolved observations are not 
possible in distant objects, and b) prove whether these fea- 
tures could be the origin —or m imic the effects — of the so- 
called temperature fluctuations (|Peimbertlll967T ). Finally, it 
must be remarked that abundance determinations based on 
recombination line ratios —which are almost independent 
the adopted n e and T e — are virtually unaffected by these 
disturbing effects. 
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